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Abstract 

The inclusive distributions of gluons and pions are calculated with 
absolute normalization for high-energy nucleon-nucleon collisions. The 
results for several unintegrated gluon distributions from the literature 
are compared. The gluon distribution proposed recently by Kharzeev 
and Levin based on the idea of gluon saturation is tested against 
DIS data from HERA. We find huge differences in both rapidity and 
transverse momentum distributions of gluons and pions in nucleon- 
nucleon collisions for different models of unintegrated gluon distri- 
butions. The approximations used recently in the literature are dis- 
cussed. The Karzeev-Levin gluon distribution gives extremely good 
description of momentum distribution of charged hadrons at midra- 
pidities. Contrary to a recent claim in the literature, we find that the 
gluonic mechanism discussed does not describe the inclusive spectra 
of charged particles in the fragmentation region, i.e. in the region of 
large \y\ for any unintegrated gluon distribution from the literature. 

1 Introduction 



The recent results from RHIC (see e.g. P) have attracted renewed interest 
in better understanding the dynamics of particle production, not only in 
nuclear collisions. 

Quite different approaches [21 El Ej have been used to describe the particle 
spectra from the nuclear collisions j^. The thermal models do not make a di- 
rect link to nucleon-nucleon collisions. In contrast, in dual parton approaches 
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(DPM) the nucleon-nucleon collisions are the basic ingredients of nuclear 
collisions. Somewhat extreme model in Ref. [3] with an educated guess for 
unintegrated gluon distribution describes surprisingly well the whole charged 
particle rapidity distribution by means of gluonic mechanisms only. Such a 
gluonic mechanism would lead to the identical production of positively and 
negatively charged hadrons. The recent results of the BRAHMS experiment 
12] put into question the successful description of Ref . and show that the 
DPM type approaches seems more correct. In the light of the BRAHMS 
experiment is becomes obvious that the large rapidity regions have more 
complicated flavour structure. The pure gluonic mechanisms, if at all, can 
be dominant only at midrapidities although the charged kaons [5J show that 
even this is doubtful. Similarly also the thermal models have difficulties to 
describe the (pseudo) rapidity dependence of particle to antiparticle ratios 
IS] and have to limit to the midrapidity only. In principle, the dynamics in 
nucleus-nucleus collision is fairly complicated and requires a separate anal- 
ysis. In the following I concentrate only on nucleon-nucleon collisions - the 
basic ingredients of the nucleus-nucleus collisions. 

On the microscopic level the approach of Kharzeev and Levin |3] is based 
on the gluon-gluon fusion. The gluon-gluon fusion is expected to be the 
dominant process at midrapidities and at asymptotically large energies. It is 
not clear how large the energy should be to validate this thesis. The physics 
in the fragmentation region is somewhat different. It was suggested long ago 
in] that pions in the fragmentation region are correlated with the valence 
quark distributions in hadrons. 

The standard hadronization approaches are based rather on the 2 
2 partonic subprocesses which constitute only a part of the dynamics. The 
perturbative component of these hybrid models has a flavour structure as dic- 
tated by the quark/antiquark distributions. On the other hand, the flavour 
structure of the remaining soft component is not so explicit. Furthermore the 
partition into "soft" and "hard" components is somewhat arbitrary, being to 
some extend rather an artifact of a natural failure in applying the (2 —>■ 2) 
pQCD at low transverse momenta of hadrons than a clear border of the two 
regions. 

In this paper I discuss the relation between unintegrated gluon distri- 
butions in hadrons and the inclusive momentum distribution of particles 
produced in hadronic collisions. The results obtained with different uninte- 
grated gluon distributions presented recently in the literature are shown and 
compared. In the present study I limit to the nucleon-nucleon collisions only 
and leave the nucleus-nucleus collisions for a separate analysis. 
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2 Photon-nucleon cross section at high ener- 
gies 

It became a standard in recent years to first describe the HERA data and 
only then to test the resulting gluon distributions in other processes. We try 
to follow this reasonable methodology also for jet and particle production. 

It is known that the LO total 7*A^ cross section can be written in the 
form 



In this paper we take the so-called quark- antiquark photon wave function 
of the perturbative form j7j. As usual, in order to correct the photon wave 
function for large dipole sizes (nonperturbative region) we introduce an ef- 
fective quark/antiquark mass {fTLeff = mo). The dipole-nucleon cross section 
can be parametrized or calculated from the unintegrated gluon distribution 



In the equation above the running coupling constant is fixed constant or is 
frozen according to an analytic prescription jSj. In the next section we shall 
compare the dipole-nucleon cross sections calculated from different uninte- 
grated gluon distributions. 

3 Unintegrated gluon distributions 

Search for the unintegrated gluon distribution in the nucleon was a subject of 
active both theoretical and phenomenological research in recent years. Still 
at present the unintegrated gluon distributions are rather poorly known. The 
main reason of the difficulties is the fact that the unintegrated gluon distri- 
bution is a quantity which depends on at least two variables [x and k^) in 
a nontrivial and a priori unknown way. Another difficulty is in an unam- 
biguous separation of perturbative and nonperturbative regions. In general 
different phenomena test the unintegrated gluon distribution in different cor- 
ners of the phase space. Therefore it is not surprising that different gluon 
distributions found in the literature, extracted from the analyses of differ- 
ent phenomena, differ among themselves considerably [0]. In this section I 
collect and briefly discuss gluon distributions used in the present calculation 
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of the jet and particle production. There are two different conventions of 
introducing unintegrated gluon distributions in the hterature. The resuhing 
quantities are denoted as / (dimenionless quantity) and T (with dimension 
1/GeV^). We shall keep this notation throughout the present paper. 



3.1 BFKL gluon distribution 



At very low x the unintegrated gluon distributions are believed to fulfil BFKL 
equation jJOl (see also PT]). After some simplifications |12. the BFKL equa- 
tion reads 



-X- 



d f{x,ql 
dx 



asN, 



dq 



fix^qft) - fix^Qt) , f{x,qt) 



qlt\ 
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(3) 



The homogeneous BFKL equation can be solved numerically J2] ■ Here in the 
practical applications we shall use a simple parametrization for the solution 

m 



J [X, K,^ ) — —r \ 
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exp 



r2A"ln(l/x) 



(4) 



In the above expression A = 4Q;5ln2, A" = 28 asC(3), cts = 3as/7r,C(3) = 
1.202. The remaining parameters were adjusted in ^3] to reproduce with 
a satisfactory accuracy the gluon distribution which was obtained in [12] as 
the numerical solution of the BFKL equation. It was found that q = qo = 1, 
C0O = 1.19 and r = 0.15 



3.2 Golec-Biernat-WiisthofF gluon distribution 

Another parametrization of gluon distribution in the proton can be obtained 
based on the Golec-Biernat-Wiisthoff parametrization of the dipole-nucleon 
cross section with parameters fitted to the HERA data [H]. The resulting 
gluon distribution reads 15]: 

a,J^(x, K^) = ^RI{x)k^ exp{-Rl{x)K^) , (5) 

where 

\ A/2 
X \ 



From their fit to the data: do = 29.12 mb, Xq = 0.41 ■ 10"^, A = 0.277 
In order to determine the gluon distribution needed in calculating jet and 
particle production we shall take a.^ = 0.2. 
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3.3 Kharzeev-Levin gluon distribution 

Another parametrization, also based on the idea of gluon saturation, was 
proposed recently in In contrast to the GBW approach [T3j, where the 
dipole-nucleon cross section is parametrized, in the Karzeev-Levin approach 
it is the gluon distribution which is parametrized. In the following we shall 
consider the most simplified functional form: 



if < Q 



The saturation momentum is parametrized exactly as in the GBW model 
Ql{x) = 1 GeV^ • (^) . It was claimed in [3j that the gluon distribution like 
dZj) leads to a good description of the recent RHIC rapidity distributions. It 
is interesting to check its performance for the deep inelastic scattering at low 
Bjorken x. 

In the following the normalization constant /o is adjusted to roughly de- 
scribe the HERA data. We find /o = 170 mb. The quality of the fit is 
shown in Fig. Q for = 0.25, 5, 10 GeV^. In this fit the running coupling 
constant frozen according to 8J was used. The result at low photon virtu- 
ality (Q^ = 0.25 GeV^) depends also on the value of the quark/antiquark 
effective mass. In the calculation in Fig. [T] mo = 0.15 GeV (sohd) and mo 
= 0.10 GeV (dashed) was used. It can be inferred from the figure that the 
(virtual) photon-proton cross section at large virtuality {Q"^ = 5,10 GeV^) is 
in practice independent of the effective quark mass. This allows to fix the 
gluon normalization constant /q. 

In order to better visualize the difference to the GBW model, in Fig. |21 1 
compare the dipole-nucleon cross sections in both parametrizations for x = 
10~^, 10~^, 10~^. In the GBW approach the dipole-nucleon cross section 
saturates at large dipole size. In contrast to the GBW parametrization, the 
dipole-nucleon cross section calculated according to Eq.Q based on the KL 
gluon distribution ((7j) grows slowly with the dipole size p. 

3.4 Kimber-Martin-Ryskin gluon distribution 

The unintegrated gluon distribution can be obtained even when the inte- 
grated gluon distribution fulfils standard DGLAP evolution equation. At 
very small x 

d 

^{x, 1^'^) = [xgix, Q^)] |q2=«2 . (8) 

This prescription breaks at larger values of x when the derivative of the 
gluon distribution becomes negative. This may be somewhat improved by 
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Figure 1: The cross section ol^f as a function of the center of mass energy 
W for g2 = 0.25 GeV^, = 5 GeV^ and = 10 GeV^ . The results 
obtained with the KL gluon distribution (niQ = 0.15/0.10 GeV) are shown by 
the solid and dashed lines. Experimental data were taken from \25^ . 
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Figure 2: T/ie dipole-nucleon cross section as a function of the transverse 
dipole size p for the GBW (dashed) and KL (dotted) unintegrated gluon dis- 
tributions. 



introducing a Sudakov form factor Tgi^n?,^'^). Then the unintegrated gluon 
distribution reads 



d 

T{x, I?) = [T(Q^ fi'^)xg{x, Q^)] |q2=^2 . (9) 

Resumming virtual contributions to DGLAP equation, the unintegrated 
parton distributions can be written as [IH] 

Ux,K',^^')=UK',^^').^^Yl r' Paa'iz){-)a'[-^,K')dz. (10) 

vr z z 

Speciahzing to the gluon distribution the Sudakov form factor reads as 



\ .L2 p^ 2tt Jq 



Tg{K^,fi'^) = exp I - / -^^v/^ J I 



(11) 

The Sudakov form factor introduces a dependence on a second scale /i^. It 
is reasonable to assume that the unintegrated gluon density given by (fTUj) 
starts only for > k^q At lower an extrapolation is needed. In 

our case of particle distributions the results are sensitive to rather low k. 
Because of this, a use of the GRV integrated gluon distribution [ini ED] in 
()10|) seems more adequate than any other PDF. Following Ref. [T7] k^q = 0.5 
GeV^ is taken as the lowest value where the unintegrated gluon distribution 
is calculated from Eg. (fTUj) . Below it is assumed 

^(x, K^) = fix, = fix, kD/kI . (12) 

The choice of /i^ in our case of jet (particle) production is not completely 
obvious. In the present analysis /i^ = pf is assumed, where pt is transverse 
momentum of the produced gluon (= jet). In accord with the interpreta- 
tion of the Sudakov form factor as a survival probability we assume that if 
transverse momentum of the produced gluon is smaller than the transverse 
momentum of the last gluon of the ladder (pt < ki or pt < 1^2, see next sec- 
tion) then the corresponding Sudakov form factor is set to 1, i.e. we do not 
allow for any enhancement. If Tg in Eq. pO|) is ignored we shall denote the 
corresponding gluon distribution as fncLAP or J-'dglap and call it DGLAP 
gluon distribution for brevity. 

3.5 Bliimlein gluon distribution 

In the approach of Bliimlein j2Tj the nf dependent gluon distribution sat- 
isfying the BFKL equation can be represented as the convolution of the 
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integrated gluon density xg^x.jj?) and a universal function B 

^{x,Ht,fi^)= B{z,K'f,fi'^)-g{-,n'^)dz . (13) 

The universal function i3(x, Kj,/i^) can be represented as a series |21j. The 
first term of the expansion describes BFKL dynamics in the double-logarithmic 
approximation: 



^^'"^'^^ \^/o(2V«.log(l/z)log(Kf//x2)) ifK?>;x2 ^ ^ 

where = Sas/ir- 

In DIS there is a natural choice of the scale /i^. The choice of the scale 
is not so obvious in the case considered in the present paper. I shall argue 
that in practice the dependence on that scale is very weak. In the following 
I shall use the integrated gluon distribution in Eq. p3|l from Ref. [TH] . 



4 Inclusive gluon production 

Before we go to particle production in the next section, let us consider the 
first step of the process - production of partons. At high energies gluons 
are the most abundantly produced partons in hadron-hadron collisions. Also 
gluons are responsible for their production. 

At sufficiently high energy the cross section for inclusive gluon production 
in /ii + /i2 ^ g can be written in terms of the unintegrated gluon distributions 
"in" both colliding hadrons: 

'J~ir~ = 1^77~T~2 / «s(^^^)^l(a;i, l4)^2ix2, i4)5{ki + K2 - Pt) d^Kid^K2 . 

dyd'^pt -^Pl J 

(15) 

In the equation above /i and /2 are unintegrated gluon distributions in 
hadron hi and /;.2, respectively. The longitudinal momentum fractions are 
fixed by kinematics: Xij2 = ^ ■ exp(±y). Generally the smaller jet (parton) 
momenta Pf, the smaller Xi/2 come into play. The argument of the running 
coupling constant is taken as VL^ = max(/t^, /t2,Pj ). The formula p5|l above 
was first written by Gribov, Levin and Ryskin (22] (see also jSHI) and used 
later e.g. in [13j. As discussed in Ref. [2^ the normalization of the cross 
section in some previous works was not always correct. Making use of the 6 
function (momentum conservation) one can simplify (fT3|) to the integral 

(16) 



dyd'^Pt -'^Pt^ 
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where Qt = k,i — ^2 was introduced. The factor 1/4 is the jacobian of trans- 
formation from (ki, R2) to {pti qt). The integral above is a two-dimensional 
integral over (Pqt, i.e. over qtdqtdcj), where is the azimuthal angle between 
qt and pf The original integral ()l(i|l can be written as 



da 



dyd^pt 

where 



rf0 , (17) 
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Figure 3: The intrinsic azimuthal correlations for different unintegrated gluon 
distributions: GBW (dashed), KL (solid), BFKL (dotted), Bliimlein (thick 
dash-dotted), DGLAP (thick dashed) and KMR (thick solid) at W = 200 
GeV. 



In Fig. ini I show the intrinsic angular correlation function /(0) for dif- 
ferent models of unintegrated gluon distributions for a RHIC energy W = 
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200 GeV. In this calculation y=0 and pt = 1 GeV was taken. Quite a dif- 
ferent pattern is obtained for different unintegrated gluon distributions. The 
0-distribution is flat for the KL, BFKL, DGLAP and KMR gluon distribu- 
tions. The most pronounced structure is obtained with the Bliimlein gluon 
distribution pT](GRV95, fi^ = 10 GeV^). It was checked that the Bliimlein 
(GRV95) gluon distribution is not very sensitive to the choice of the second 
scale fj,"^. The dependence ai y ^ also strongly depends on the uninte- 
grated gluon distribution. 

It was suggested in |3] that the integral (fTHjl may be approximated by the 
formula 

^ ' (19) 




Figure 4: A comparison of the gluon rapidity distributions obtained from 
the exact / Ii6|) (thick lines) and approximate My\) (thin lines) formula for 
different models of unintegrated gluon distributions at W = 200 GeV. 



In Figl^ (rapidity distribution) and FigEl (transverse momentum distri- 
bution) I compare the results using the exact Eq. p(j|l and the approximate 
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Figure 5: A comparison of the gluon transverse momentum distributions 
obtained from the exact M(^) (thick lines) and approximate M9) (thin lines) 
formula for different models of unintegrated gluon distributions at W = 200 
GeV. 



Eq.()19p formulae for different models of unintegrated gluon distributions. In 
FigEJthe integration over pt > 0.5 GeV is performed while in Fig|Hl-l < y < 
1. In both cases as was fixed at 0.2. As can be seen by inspection of the 
figures the use of the approximate formula is quantitatively justified for the 
KL, BFKL gluon distributions and not justified for the GBW one. 

In Figinil compare the cross section ^{y) for different models of uninte- 
grated gluon distributions. In this calculation pt < 0.5 was assumed. The 
rapidity distribution of gluons are rather different for different gluon PDF. 
Average values of xi and X2 obtained with different gluon distribution with 
the Pt interval chosen are shown in FigI7| The following general observations 
can be made. Average value < Xi > and < X2 > only weakly depend on 
the model of unintegrated gluon distribution. For ?/ ~ at the RHIC en- 
ergy W = 200 GeV one tests unintegrated gluon distributions at Xg = 10~^ 
- 10^^. This is the region known already from the HERA kinematics. When 
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Figure 6: Inclusive gluon rapidity distribution (pt > 0.5 GeV) at W = 200 
GeV for different models of unintegrated gluon distributions. 
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\y\ grows one tests more and more asymmetric (in xi and X2) configurations. 
For large \y\ either xi is extremely small {xi < 10"^) and X2 — * 1 or a;i — > 1 
and X2 is extremely small (x2 < 10~^). These are regions of gluon momen- 
tum fraction where the unintegrated gluon PDF is rather poorly known. The 
approximation used in obtaining unintegrated gluon distributions are valid 
certainly only for x < 0.1. In order to extrapolate the gluon distribution 
to — 1 I multiply the gluon distributions from the previous section by a 
factor (1 — Xg)^, where n = 5-7. 




-5 -2.5 2.5 5 

y 

Figure 7: The average value of Xi and X2 forpt > 0.5 GeV and at W = 200 
GeV. Lines corresponding to different unintegrated gluon PDF are identical 
as in the previous figure. 

In the approach considered in the present paper (for details see next sec- 
tion) the production of particles is sensitive to rather small gluon (called 
equivalently jet despite of the small transverse momentum) transverse mo- 
menta. 

In FigUl plot ^(pt) in the low Pt region. In these calculations the gluon 
rapidity was integrated in the interval -1 < y < 1. The results obtained with 
different models for unintegrated gluon distributions differ considerably. The 
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Figure 8: Inclusive gluon transverse momentum distribution (-1 < y < 1) 
at W = 200 GeV for different models of unintegrated gluon distributions: 
BFKL (dotted), GBW (dashed), KL (solid), Bliimlein (thick dash- dotted) 
and DGLAP (thick dashed). 



transverse momentum distribution obtained with the GBW gluon density is 
much steeper than the distribution for any other gluon density. The inclu- 
sion of DGLAP evolution as in [211 would probably change the situation. In 
the case of the Bliimlein gluon distribution the transverse momentum spec- 
trum has a natural low-pt cut-off if the scale /i^ = pi is chosen. If similar 
prescription of the scale is used for calculating gluon transverse momentum 
distribution with KMR method the DGLAP and KMR results are almost 
identical. Contrary to the claim in P] the result obtained with the GBW 
and KL gluon distributions differ considerably. 

The rapidity and pseudorapidity distributions of partons (massless par- 
ticles) are identical. The situation changes when massive particles are pro- 
duced in the final state via fragmentation. Below we discuss how to take into 
account the unknown hadronization process with the help of phenomenolog- 
ical fragmentation functions. 
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5 From gluon to particle distributions 

In Ref. jS] it was assumed, based on the concept of local parton-hadron du- 
ality, that the rapidity distribution of particles is identical to the rapidity 
distribution of gluons. This seems to be a very severe assumption and for 
massive particles this idea must lead to incorrect results, especially in the 
fragmentation region. This approach leads to e.g. (massive) particles with 
rapidities (i/h) beyond the allowed kinematical region {yh,min,yh,max)- Fur- 
thermore in [3] the normalization of rapidity distributions was fitted to the 
experimental charged particle rapidity distributions. In our opinion, the good 
description of the charged particle distribution in the full range of rapidity 
in Ref. ^ is due to these simplifications rather than due to the underlying 
dynamics. 

In the present approach I follow a different, yet simple, approach which 
makes use of phenomenological fragmentation functions (see e.g.|27^ ES!)- 
For our present exploratory study it seems sufficient to assume that the 
emitted hadron, mostly pion, is collinear to the gluon direction {6^ = Og). 
This is equivalent to f]h = Vg = Vgy where rjh and rjg are hadron and gluon 
pseudorapitity, respectively. 

In experiments a good identification of particles is not always achieved 
which makes impossible to determine the rapidity of a particle. The practice 
then is to measure pseudorapidity. The rapidity of a given type of hadrons 
{yh) with a mass nih can be obtained from the pseudorapidity as 



The collinearity of partons and particles leads to the following relation be- 
tween rapidity of the gluon and hadron 



nomenological fragmentation functions one has to define a new kinematical 
variable. In accord with e~^e~ and ep collisions I define a standard auxiliary 
quantity z by the equation Eh = zEg. This leads to the following relation 
between transverse momenta of the gluon and hadron 




(20) 




(21) 




In order to introduce phe- 



Pt,g 



Pt,h 



J{mt,h,yh) , 



(22) 



z 
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where 

J{mt,H,yH)=(l- , ) . (23) 

\^ mf^^ cosh yh J 

Now we can write the single particle distribution in terms of the gluon dis- 
tribution from the last section as follows 

da{r]h,pt,h) 



dr]hd'^Pt,h 
KVg - Vh) 5^ ( Vt,h - 



dygd pt^g J dz Dg_^h{z, /i£,) (24) 
zpt,g\ da{yg,pt,g) 



J J dygd^pt^g 

Making use of the 5 functions we can write the single particle spectrum as 



dai'nh, Pt,h) _ f'"'"^ ^^J^Dg_h{z, /il)) da{yg, pt^. 



dVhd'^Pt,h Jz^,^ z"^ dygd^pt^g 



(25) 



Experimentally instead of the two-dimensional spectrum ()25|) one determines 
rather one-dimensional spectra in either rj^ or pt^^. 

The one-dimensional pseudorapidity distribution can be obtained by in- 
tegration over hadron transverse momenta 

d(y{rih) f ,2 da{rih,pt,h) 

^ " Pt,h —3 — 15 • (26) 



dVh J ' drjhd^pt^h 

Stable particles ^ are produced directly in the fragmentation process or 
are decay products of other unstable particles. There are a few global anal- 
yses of fragmentation function in the literature up to next-to-leading order 
IHOl inn inn]- in the present calculation I shall use only leading order 
fragmentation functions from |3ni |M] . One should remember, however, that 
both e"'"e~ and ep collisions do not allow to uniquely determine Dg^h frag- 
mentation functions. In order to test sensitivity of our results to these, in 
my opinion, not quite well known objects I shall use also simple functional 
forms: Dg^h{z) =2^ (model I) or Dg^h{z) = 3^^^ (model II) with the 
factors in front adjusted to conserve momentum sum rule. When charged 
particles are measured only, then to a good approximation it is sufficient to 
multiply the fragmentation functions above by a factor 2/3. 

In Figini I compare pseudorapidity distribution of charged pions at W = 
200 GeV calculated with the KL gluon distribution and different parametriza- 
tions of fragmentation functions. For the BKK1995 30j and for the KKP2000 

^Here by stable particles we mean the particles registered in detectors 
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Figure 9: Charged-pion pseudorapidity distribution at W = 200 GeV for the 
KL unintegrated gluon distribution for different parametrizations of fragmen- 
tation functions. In this calculation pt^h > 0.2 GeV. The experimental data 
of the UA5 collaboration are taken from IS^ . 



\M\ fragmentation functions the factorization scale was set to ~ Ptg^ 
cept for g < 1 GeV, where it was frozen at = 1 GeV^. For reference 
shown are also experimental data for charged particles measured by the UA5 
collaboration at CERN [33] . The results only weakly depend on the choice of 
the g ^ TT fragmentation function. It is worth stressing that the theoretical 
cross section at rjh ^ is almost consistent with the experimental one. How- 
ever, the shapes of theoretical and experimental pseudorapidity distributions 
differ significantly. It seems there is a room for different mechanisms typical 
for fragmentation regions. The specificity of these regions will be discussed 
elsewhere. 

Let us analyze now how the results for pseudorapidity distributions de- 
pend on the choice of the unintegrated gluon distribution. In FigEH I com- 
pare pseudorapidity distribution of charged pions for different models of 
unintegrated gluon distributions. In this calculation the Binnewies-Kniehl- 
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Figure 10: Charged-pion pseudrapidity distribution at W = 200 GeV for 
different models of unintegrated gluon distributions. In this calculation pt,h > 
0.2 GeV. The experimental data of the UA5 collaboration are taken from Vdd^ . 



Kramer fragmentation function [20] has been used. The conclusions inferred 
above stay true also here. Having in view a dramatically steep pt^g distribu- 
tion in FiglHl it is rather surprising that the normalization of the spectra at 
midrapidities comes roughly correct, although very is a tendency to an over- 
estimation for some gluon distributions. This can be due to the fit to DIS 
data, where the resolved photon component has been neglected. If the re- 
solved photon component is explicitly included jHl] then the normalization of 
the dipole-nucleon component (dipole-nucleon cross section or unintegrated 
gluon distribution) must be reduced. 

What are typical transverse momenta of gluons involved in the calcula- 
tions is shown in Fig|TT] In this calculation we have used the KL unintegrated 
gluon distribution and the BKK g ^ ir fragmentation functions [HOI- We ob- 
serve a maximum of the transverse momentum squared of the produced gluon 
at Tjh ~ 0. In our implementation of fragmentation {pjj^ <^ p^^g) one tests 
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Figure 11: Average values of {pf) (solid), {k{) (dashed) and (dotted) as 
a function of pion pseudorapidity 



relatively large pfg. While at midrapidities {pf) > {i^i)^ {1^2)1 when going 
to the fragmentation regions the relation reverses. In the whole range of 
pseudorapidity one tests on average ~ 1 GeV^. One should remem- 

ber, however, that at the same time (xi) and (^2) change dramatically when 
going from midrapidities to the fragmentation region. 

In contrast to Ref.|2j, where the whole pseudorapidity distribution, in- 
cluding fragmentation regions, has been well described in an approach simi- 
lar to the one presented here, in the present paper pions produced from the 
fragmentation of gluons in the gg ^ g mechanism populate only midrapid- 
ity region, leaving room for other mechanisms in the fragmentation regions. 
These mechanisms involve quark/ ant iquark degrees of freedom or leading 
protons among others. In Fig|T21 1 show the pseudorapidity spectra of pro- 
tons, antiprotons and the difference da/driT^+ — da/drj^^- obtained with the 
code HIJING (see also [38]). The difference of the proton-antiproton 
spectra gives an idea of leading particle contribution. Both protons from 
deeply inelastic events as well as protons from diffraction dissociation (sin- 
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Figure 12: The pseudorapidity distribution of protons (solid), antiprotons 
(dashed) and the difference of the spectra of and t^^ (dash-dotted) in the 
proton-proton collision at W = 200 GeV obtained with the code HIJING 
The thick solid line corresponds to the sum of these three contributions. The 
experimental data of the UA5 collaboration are taken from 



gle diffraction) have been included. The difference of the positively and 
negatively charged pions gives the lower limit on the vr"*" — vr" asymmetric 
mechanisms not taken into account in the Kharzeev-Levin approach. The 
sum of the three contributions (thick solid) gives then lower limit on the 
missing contributions. It is of the similar size as the missing contributions 
in FigEland FiglTUl This strongly suggests that the agreement of the result 
of the gg ^ g approach with the PHOBOS distributions j3] in Ref.j2] in the 
true fragmentation region is rather due to approximations made in |3j than 
due to correctness of the reaction mechanism. In principle, this can be veri- 
fied experimentally at RHIC by measuring the tt^/tt^ ratio in proton-proton 
scattering as a function of (pseudo) rapidity in possibly broad range. It seems 
that the BRAHMS experiment, for instance, can do it even with the existing 
apparatus. 
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Figure 13: Transverse momentum distributions of charged pions at W = 200 
GeV for the KL gluon distribution and different fragmentation functions. 
The experimental data of the UAl collaboration are taken from 15'6) /. 



The transverse momentum distribution of charged hadrons is shown in 
FigEl together with experimental data of the UAl collaboration at CERN 
from Ref ■ [SB] . In this calculation the KL gluon distribution has been used. 
It is not completely clear to me how the experimental data in jSQl should be 
interpreted. ^ I assume that the experimental data should be interpreted as: 

X = [ j^dv, / [dv,. (27) 
J drjhd^pt J 

We have taken rj^ G (-2.5,2.5). The simple hadronization functions, called 
model I and II above, correctly fit low pt^h data and fail in the large pt^h 
region. This is due to lack of QCD evolution jSHj- The results obtained 
with fragmentation functions from l!TT] which include DGLAP evolution, 
extremely well describe the large pt^h data. Having in mind the ambiguity of 

■^The notion of the invariant cross section in |36j is contradictory to the lack of particle 
identification there. 



22 



the experimental data interpretation, the KL gluon distribution does a fairly 
good job. 

In Fig^^ I compare the theoretical transverse momentum distributions 
of charged pions obtained with different gluon distributions with the UAl 
collaboration data (^EI- The best agreement is obtained with the Karzeev- 
Levin gluon distribution. The distribution with the GBW model is much too 
steep in comparison to experimental data. This is probably due to neglecting 
QCD evolution. 




Figure 14: Transverse momentum distributions of charged pions at W = 200 
GeV for BKK1995 fragmentation function and different models of uninte- 
grated gluon distributions. The experimental data of the UAl collaboration 
are taken from 



6 Conclusions 

I have calculated the inclusive distributions of gluons and associated charged 
pions in the nucleon-nucleon collisions through the gg ^ g mechanism in the 
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/cf-factorization approach. The results for several unintegrated gluon distri- 
butions proposed recently in the literature have been compared. The results, 
especially transverse momentum distributions, obtained with different mod- 
els of unintegrated gluon distributions differ considerably. 

A special attention has been devoted to the gluon distribution proposed 
recently by Kharzeev and Levin to describe charged particle production in 
relativistic heavy-ion collisions. In the first step I have tested the gluon dis- 
tribution in electron deep-inelastic scattering at small Bjorken x. A rather 
good description of the HERA data can be obtained by adjusting a normal- 
ization constant. In the next step so-fixed gluon distribution has been used to 
calculate (pseudo)rapidity and transverse momentum distribution of gluonic 
jets and charged particles. 

Huge differences in both rapidity and transverse momentum distributions 
of gluons and pions for different models of unintegrated gluon distributions 
have been found. 

Some approximations used recently in the literature have been discussed. 
Contrary to a recent claim in Ref . [3] , we have found that the gluonic mech- 
anism discussed does not describe the inclusive spectra of charged particles 
in the fragmentation region, i.e. in the region of large (pseudo)rapidities for 
any unintegrated gluon distribution from the literature. Clearly the gluonic 
mechanism is not the only one and other mechanisms (see e.g.pTf EHj ) ne- 
glected in [5] must be added. Some of them have been estimated with the 
help of the HIJING code, giving a right order of magnitude for the missing 
strength. 

Since the mechanism considered is not complete, it is not possible at 
present to precisely verify different models of unintegrated gluon distribu- 
tions. The existing gluon distributions lead to the contributions which almost 
exhaust the strength at midrapidities and leave room for other mechanisms 
in the fragmentation regions. It seems that a measurement of transverse mo- 
mentum distributions of particles at RHIC should be helpful to test better 
different unintegrated gluon distributions. A good identification of particles 
is required to verify the other mechanisms. 

In contrast to standard integrated gluon distributions, the extraction of 
unintegrated gluon distribution from experimental data seems a rather diffi- 
cult task. At present, one can rather test different unintegrated gluon distri- 
butions based on different models existing in the literature. In the present 
analysis I have discussed whether the production of particles can provide 
some information on unintegrated gluon distributions in the nucleon. There 
are many other reactions where this is possible, to mention here only heavy 
quark or jet production in ep and pp collisions. Going to more exclusive 
measurements seems indispensable. An example is a careful study of jet 
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correlation in photon-proton ^39j and nucleon-nucleon [AQ^ collisions. In my 
opinion, we are at the beginning of the long way to extract gluon or more 
generally parton unintegrated distributions. 
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